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Abstract. Despite the importance of Type la supernovae as standard candles for cosmology and 
to the chemical evolution of the Universe, we still have no consistent picture of the nature of these 
events. Much progress has been made in the hydrodynamical explosion modelling of supernovae 
la in the last few years and fully 3-D explosion models are now available. However those simula- 
tions are not directly comparable to observations: to constrain explosion models, radiative transfer 
calculations must be carried out. We present a new 3-D Monte Carlo radiative transfer code which 
allows forward modelling of the spectral evolution of Type la supernovae from first principles, using 
hydrodynamical explosion models as input. Here, as a first application, we calculate line-of-sight 
dependent colour light curves for a toy model of an off-centre explosion. 
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INTRODUCTION 

The greatest predictive power of explosion models for Type la supernovae can be 
extracted by performing parameter free radiative transfer calculations to directly link 
them to observational data. This requires a solution of the time-dependent 3-D radiative 
transfer problem in chemically inhomogeneous models of supernova ejecta. 

METHOD 

We extended the 3-D Monte Carlo radiative transfer code introduced in [Q]] to a non^grey 
opacity treatment following the ideas outlined in a series of papers by Lucy [0, SH. 
The basic assumption of this approach is that the ejecta are in homologous expansion 
thereby decoupling the radiative transfer from hydrodynamical explosion modelling. 
As input we therefore take densities, velocities and composition from explosion models 
which have been extended up to the phase of homologous expansion and map these 
on to a 3-D cartesian grid which expands with time. The total energy released by the 
synthesised 56 Ni is divided into N identical energy packets which are distributed on the 
grid according to the 56 Ni distribution. These pellets follow the homologous expansion 
until they decay. Decay times are sampled randomly according to the 56 Ni^ 56 Co— > 56 Fe 
decay chain. 

Upon decay, a pellet transforms to a y-packet representing a bundle of monochromatic 
y-radiation. The y-packets then propagate through the grid until they either interact 
with matter by Compton scattering, photoabsorption or pair creation or they leave the 



simulation volume. When a y-packet interacts it can transfer its energy to an electron. 
We assume that the timescale to thermalise these fast electrons is short and convert 
the y-packet into a packet of thermal kinetic energy, a ^-packet in the nomenclature of 
Lucy QBHl- ^-packets are not propagated but are converted in situ to radiative energy 
by sampling all the possible cooling channels, i.e. collisional excitation/ionisation, 
bound-free and free-free cooling. 

In the case of bound-free or free-free cooling the ^-packet is transformed into a 
monochromatic energy packet representing ultraviolet-optical-infrared (UVOIR) ra- 
diation, a so-called r-packet. The frequency of the packet is sampled randomly in 
accordance with the selected cooling process, r-packets propagate through the grid 
either until they interact with matter by electron scattering, free-free absorption, bound- 
free absorption or line absorption or until they leave the simulation volume. For free-free 
absorptions the r-packets transform into fc-packets, for line absorptions they transform 
into a packet of atomic internal energy while for bound-free absorption they may be 
converted to either type. Atomic internal energy is instantaneously transformed into 
thermal kinetic energy by collisional deexcitation/recombination or radiative energy by 
line emission/radiative recombination by sampling the statistical equilibrium equations 
following the macro-atom approach of Lucy Q2Q . 

When the simulation has finished, we extract the spectral evolution by binning the 
escaping r-packets by frequency, time and angle. Colour light curves are extracted from 
the spectral evolution by integrating the spectra over the filter functions. 

In calculating opacities we need to know atomic level populations and therefore the 
radiation field. In principle, these could be extracted exactly out of our simulation, but it 
is too computationally demanding to do so. Therefore, we parameterise the radiation 
field in all grid cells in the nebular approximation using radiation temperatures and 
dilution factors from Monte Carlo estimators (c.f. Mazzali & Lucy [50). The ionisation 
balance is solved exactly using Monte Carlo estimators for the ionisation rates and 
the excitation state is computed assuming the Boltzmann distribution evaluated at a 
temperature corresponding to the local energy density of the radiation field. Full details 
will be presented in a forthcoming paper. 

APPLICATION TO AN OFF-CENTRE TOY MODEL 

After testing the code against other radiative transfer codes [@, 0] we have applied it to 
a parameterised off-centre explosion model (like those in [|8|]) to make a preliminary 
study of line-of- sight effects on colour light curves. Our toy model has a total mass of 
1 .4 M and uniform mass density. All the matter is confined in a ball with maximum 
velocity 1.2 ■ 10 4 kms" 1 . The 56 Ni (0.52 M Q in our model) is located in an inner ball. 
This Ni-rich ball is offset from the centre of mass of the model by 2.4 • 10 3 kms -1 . Its 
extension is chosen such that the 56 Ni mass fraction inside the inner ball is 0.8. 

We mapped this model onto a 50 3 grid and followed the evolution of 10 7 packets 
over 50 timesteps from 4 to 60 days after explosion. To reduce computational costs, we 
calculated the level populations in LTE and parameterised the radiation field by a black 
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FIGURE 1. UVOIR bolometric (top left panel), U (middle left), B (middle right), V (bottom left) and R 
(bottom right) light curves for the off-centre toy model. Dashed lines show light curves as seen from the 
Ni-rich side, solid lines those seen from the Ni-poor side. The top right panel shows the 56 Ni distribution 
in the model's x-z plane. The black area contains 56 Ni, the grey has no 56 Ni. The arrows indicate the two 
lines-of-sight shown in the light curve plots. 



TABLE 1. Peak parameters of the colour light curves. 







U 


B 


V 


R 


Ni-rich side: 


Peak magnitude 


-19.58 


-18.94 


-19.16 


-18.97 




Peak time (d) 


15.30 


19.24 


19.24 


19.24 


Ni-poor side: 


Peak magnitude 


-19.27 


-18.76 


-19.01 


-18.80 




Peak time (d) 


16.02 


19.24 


21.09 


20.14 



body adopting excitation and radiation temperatures corresponding to the local energy 
density of the radiation field, which is extracted from a Monte Carlo estimator. We used 
atomic data from Kurucz CD 23 [9]. Approximately ~ 3 • 10 5 atomic lines were used. 

U, B, V, R and £/V0/i?-bolometric light curves for this model are shown in Figure 1 
along the two most extreme lines-of-sight, i.e. from the side in which the Ni bubble is 
displaced and from the opposite side. The light curves for intermediate lines-of-sight 
vary smoothly between those extreme cases. As one expects the supernova appears 
brighter on the Ni-rich side. The effect varies between 0.3 mag in [/-band to 0.15 
mag in V-band. Furthermore light curves peak earlier on the Ni-rich side (see Table 
1 for details). This confirms the trend seen in grey calculations and Sim et al. in 
this volume) even though the effects are somewhat weaker in our calculations. Further 
studies are necessary to probe a broader parameter range and more realistic models as 
well as more complete atomic data sets. 

CONCLUSION 

We presented a new time-dependent 3-D Monte Carlo radiative transfer code which 
allows forward modelling of the spectral evolution of Type la supernovae from first 
principles, using hydrodynamical explosion models as input. As a first application we 
showed line-of- sight dependent effects on colour light curves for a toy model of an off- 
centre explosion. In future work we will extend this study to more realistic models and 
the testing of state-of-the-art explosion models. 
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